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ABSTRACT 

”  A  new  model  of  the  height  profile  of  tropospheric 

refractivity  is  presented  as  the  basis  for  an  improved  tropospheric 

• 

correction  for  satellite  doppler  or  range  data.  The  model  treats 

/  1 

•;  the  "dry"  and  "wet"  components  of  the  refractivity  separately,  each 
as  a  fourth- degree  function  of  height  above  the  geoid.  The  height 
parameters  are  different  for  the  two  and  are  obtained  from  meteor¬ 
ological  balloon  data.  A  latitude  dependence  has  been  found  for  the 

dry  height.  The  resulting  expressions  for  the  correction  of 
• 

satellite  doppler  and  range  data  are  presented,  and  their  effective¬ 
ness  is  shown  by  figures  that  give  two  different  kinds  of  observed 
data:  doppler  residuals  for  several  satellite  passes,  and  the 

r  • 

navigation  error  in  station- to- orbit  slant  range.  In  each  case  the 
data  are  given  without  and  with  the  correction.  The  use  of  the 
correction  removed  obvious  systematic  errors. 
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ABSTRACT 

A  new  model  of  the  height  profile  of  tropospheric  refractivity  N 
is  presented  as  the  basis  for  an  improved  tropospheric  correction  for 
satellite  doppler  or  range  data.  (N  =  10^  (n-l),  where  n  is  the  index 
of  refraction.)  The  model  treats  the  "dry"  and  "wet"  components  of  N 
separately,  and  represents  each  as  a  fourth-degree  function  of  height 
above  the  geoid;  each  component  profile  starts  with  its  locally  observed 
surface  value  and  decreases  to  zero  at  an  effective  height  which  is  differ¬ 
ent  for  the  two  components .  The  height  parameters  were  obtained  by  a 
least  squares  fit  to  observed  data.  A  latitude  dependence  has  been  found 
for  the  "dry"  height .  The  model  was  found  capable  of  closely  matching 
any  local  average  N  profile  observed  in  a  world -wide  sample  of  locations, 
throughout  the  height  range  of  meteorological  balloon  data  (up  to  2k  km); 
samples  are  shown.  The  resulting  expressions  for  the  correction  of  satel¬ 
lite  doppler  and  range  data  are  presented;  they  are  finite  and  usable  at 
all  elevation  angles .  The  effectiveness  of  the  correction  is  evidenced  by 
figures  showing  two  different  kinds  of  observed  data:  first,  doppler 
residuals  for  several  satellite  passes,  respectively  without  and  with  the 
use  of  the  correction;  and  second,  the  "navigation"  error  in  station -to- 
orbit  slant  range  from  doppler  data,  again  without  and  with  the  correction. 
The  use  of  the  correction  removed  obvious  systematic  errors.  The  fact  that 
satellite  doppler  data  display  identifiable  tropospheric  effects  is  of 
interest  with  regard  to  future  study  of  the  troposphere. 
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I.  Introduction 

Dm  dopplcr  shift  of  a  radio  signal  frcm  a  satellite,  aa  received 
at  known  tracking  station  positions,  Is  used  In  the  Navy  navigational 
satellite  eye  tea  to  determine  the  ort>lt  of  the  satellite;  the  position  of 
any  other  observing  station  can  then  be  determined  when  the  orbit  Is  known 
[Ouler  and  USlffenbach,  I960]. 

The  apparent  (electromagnetic)  range  to  the  satellite  from  a 
tracking  station  Is  Jnds  along  the  signal  path;  n  being  the  Index  of 
refraction  of  the  medium,  varying  with  position  along  the  path.  The 
observed  dopplcr  shift  of  the  satellite  signal,  neglecting  non-atmoepherlc 
errors,  la  then 

^ob.  *  'I  It  t***  M 

where  f  Is  the  signal  frequency  and  c  the  vacuum  velocity  of  light. 

But  the  "vacuum"  dopplcr  shift  Is  needed  for  computing  the  orbit: 


“v  ■  -I  If  <2> 

where  p  Is  the  geometric  slant  range  to  the  satellite. 

Both  Ionospheric  and  tropospheric  effects  must  be  removed  from 
the  observed  doppler  to  get  the  vacuum  dopplcr  shift.  The  two-frequency 
method  [Ouler  and  Weiffenbach,  i960]  which  removes  first-order  ionospheric 
effects  does  not  remove  the  effects  of  the  troposphere,  whose  refractive 
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Index  Is  independent  of  frequency  up  to  approximately  15000  Mhz .  A 
theoretical  tropospheric  correction  is  therefore  required  and  must  be  baaed 
on  some  model  of  the  troposphere;  the  term  "troposphere"  will  be  used  here 
to  include  all  the  lower,  uncharged  part  of  the  atmosphere. 

An  initial  tropospheric  mcctel  and  a  subsequent  improvement 
[Hopfleld,  1963,  1965], which  will  here  be  designated  as  Model  1  and 
Model  II,  have  been  euperseded  by  the  model  to  be  deacrlbed  below. 

n.  Assumptions  and  Requirements  for  a  Troposphere  Model 

It  will,  be  assumed,  as  before,  that  the  refractive  index  of  the 
tropoephere  is  a  function  of  its  value  at  the  surface  and  of  height  above 
the  earth,  but  not  of  horizontal  distance  or  time,  within  the  time  and 
space  pertinent  -o  a  satellite  pass.  The  curvature  of  the  signal  path 
affects  doppler  data  only  if  it  is  great  enough  t-  affect  the  path  length 
perceptibly,  therefore  only  near  the  horizon.  If  this  effect  is  neglected, 
the  Integral  of  equation  (l)  can  be  taken  along  the  slant  range  vector,  like 
that  of  equation  (2).  For  ionosphere -corrected  data  containing  only  atmos¬ 
pheric  errors,  the  difference  between  these  two  Is  the  tropospheric  effect. 
This  difference,  expressed  In  terms  of  the  refractlvlty  If  ,  is 

*tro  •  '  10-6  I  £  J*  •»  <3) 

Effectively,  the  Integral  now  extends  only  through  the  troposphere,  since 
N  becomes  zero  at  the  "top"  of  the  troposphere.  Expresolon  (3)  provides 
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the  needed  tropospheric  dppplsr  correction,  given  an  adequate  expression  for 
N  as  a  function  of  the  height  h  j  l.e.,  an  adequate  H  profile. 

To  give  viie  proper  correction  throughout  a  satellite  pass,  a 
theoretical  M  profile  aust  yield  the  observed  value  of  jfedp  through 
the  troposphere  at  all  elovatlon  angles  which  are  to  be  used.  To  do  this 
at  high  angles,  the  senlth  value,  Jfcdh  ,  aust  be  correct.  To  succeed  at 
low  angles,  tbs  theoretical  'I  profile  aust  also  natch  tbs  observed  profile 
In  shape  near  the  surface  of  the  earth.  Since  low  angle  satellite  data  can 
often  be  discarded,  tbs  fornsr  at  these  two  conditions  (tbs  zenith  Integral) 

Is  tbs  aore  essential  one  for  satellite  purposes,*  and  Is  the  condition  on 
which  Models  I  and  XX  of  the  correction  were  based.  Both  of  those  versions 
were  expected  to  be  less  satisfactory  at  low  angles  (e.g.,  below  10°)  than 
at  high  angles.  However,  both  versions  nude  equation  (3)  lntegreble  la 
closed  fora  and  both  could  be  used  even  at  tbs  horizon,  though  with  seas 
loss  of  accuracy. 

To  obtain  an  observed  N  profile  for  coaparlson  purposes,  tbs 
values  of  I  at  different  heights  were  coaputed  froa  a  set  of  asteoro- 
logical  balloon  data  [Monthly  Cllaatle  Data  tor  tbs  World],  using  tbs 
equation  of  Smith  and  Welntraub  [1993]  at  each  observed  data  point; 

*  ■  00 

K  K 

In  this,  Tg  Is  temperature  (degrees  Kslvln) ,  P  Is  total  atmospheric 
pressure  and  e  the  partial  pressure  of  water  vapor  (both  In  millibars) . 

* 

A  slallar  condition  for  the  problea  of  astroncsdcal  refraction  was  recognized 
almost  two  centuries  ego  by  Orlanl  [Mahan,  1962] 
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The  area  under  the  resulting  N  vs  h  profile  then  provides  the  observed 
value  of  the  zenith  integral  jNdh  for  the  set  of  observed  data.  The 
theoretical  profile  to  be  coapared  with  this  is  based  on  the  sane  surface 
value  of  N  ,  l.e.,  the  one  provided  by  the  observed  data. 

Model  11  of  the  theoretical  correction  succeeded  In  mateMng  a 
world -vide  sampling  of  observed  values  of  J^dh  generally  with  an  accuracy 
of  2  or  3  percent  (using  a  monthly  average  N  profile  In  each  case).  This 
version,  therefore,  '/as  expected  to  be  satisfactory  except  at  "low*'  angles, 
the  definition  of  "low"  being  sonevhat  uncertain. 

To  clarify  this  point  and  to  evaluate  the  correction,  sene 
station  positions  were  determined  from  satellite  doppJer  data,  respectively 
without  and  with  the  use  of  the  Model  II  tropospheric  correction.  Re¬ 
fraction  by  the  troposphere  steepens  the  slope  of  the  doppler  shift  vs.  time 
curve  as  observed  during  a  satellite  pass.  If  no  horizontal  gradient  of 
R  is  present  and  if  tracking  data  are  symmetrical  about  the  point  of 
closest  approach,  the  troposphere  affects  only  the  slant  range  determination 
(tracking  station  to  orbit) ,  maklr  3  the  station  appear  closer  to  the  orbit 
than  It  actually  is.  This  effect  Is  a  function  of  pass  geometry;  It  averages 
approximately  20  meters  for  a  high -angle  pass,  but  100  meters  for  a  pass 
whose  maximum  elevation  Is  only  10°.  (An  Illustration  of  the  effect  will 
be  shown  later.)  The  use  of  the  above  correction  removed  this  effect  from 
high -angle  passes,  on  the  average,  as  far  as  could  be  detected;  but  sys¬ 
tematically  over -corrected  by  approximately  10  percent  in  low-angle  passes. 
When  used  for  a  10°  puss,  it  made  the  average  station  appear  10  or  15  meters 
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too  far  away  from  the  orbit.  This  low  angle  Inaccuracy  called  for  an  Improve* 
ment  In  the  N  profile  shape. 

Ill,  Matching  the  Shape  of  the  Observed  N  Profile:  a  Two -Quart  lc 
Expression 

The  desired  theoretical  N  profile  would  start  with  the  observed 
N  at  the  surface  and  yield  not  only  the  correct  zenith  integral  jNdh  but 
also,  If  possible,  the  correct  lapse  rate  near  the  surface  of  the  earth, 
at  any  geographical  location.  Some  expressions  which  were  tried  and  dis¬ 
carded  will  be  mentioned  briefly. 

The  well-known  CRPL  exponential  atmosphere  [Bean  and  Thayer,  1959] 
gives,  to  a  good  approximation,  the  observed  lapse  rate  of  N  near  the 
surface .  It  was  designed  to  do  this  especially  in  the  continental  U.  S .  and 
was  Intended  for  low-angle  applications,  especially  within  the  troposphere . 

It  was  not  designed  specifically  for  satellite  applications  and  Its  zenith 
Integral  does  not  match  the  observed  one  for  all  climates.  If  applied  to 
marine  or  other  humid  climates,  It  Is  likely  to  make  this  Integral  too 
small,  often  by  10  or  15  percent.  In  addition,  it  Is  inconvenient  for  use 
with  doppler  data.  Since  the  integral  jNdp  for  exponential  N  ,  except 
at  the  zenith.  Is  not  integrable  In  closed  form  and  requires  approximation 
[Freeman,  1962],  the  evaluation  of  its  time  derivative  for  equation  (3)  also 
entails  approximations;  their  validity  Is  a  function  of  the  elevation  angle. 

A  model  more  suited  to  satellite  data  was  therefore  desired. 


In  contrast,  any  N  profile  of  the  form 
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N  -  k  (ho  -  h)^  ,  hlhc 

M.  t* 


(5) 


will  make  equation  (3)  integrable;  p  here  Is  a  positive  Integer  and  hQ 
the  corresponding  equivalent  height  of  the  troposphere.  A  height  h 

°M> 

can  easily  be  chosen  for  any  value  of  p  ,  to  yield  the  correct  value  of 
Jtodh  ,  if  a  correction  factor  is  used  as  for  Model  II;  but  such  a  profile 
will  not  in  general  match  the  observed  N  profile  near  the  earth.  A  simple 
profile  of  the  form  of  equation  (5),  therefore,  was  also  ruled  out. 

Separate  treatment  of  the  so-called  "dry  term"  and  "wet  term" 
of  the  refractlvity  was  then  considered.  These  are  the  components  of  N 
in  equation  (k): 


M  -  77-6  p 

‘  TK 


N 


77.6  (U810  e) 
v  T  2 
aK 


\ 


(6) 


where  d  and  w  refer  to  "dry"  and  "wet"  respectively.  Even  a  cursory 
look  at  upper  air  meteorological  data  [Monthly  Climatic  Data  for  the  World] 
shows  that  at  heights  of  9  or  10  km,  where  atmospheric  pressure  is  still 
at  nearly  one -third  of  its  surface  value,  the  moisture  content  is  nearly 
zero  (barely,  if  at  all,  measurable),  even  in  tropical  latitudes.  The 
components  Nd  and  Nw  are  therefore  different  functions  of  height. 

Evidence  that  only  a  part  of  the  tropospheric  range  or  range  rate 
error  is  weather  dependent  had  been  found  in  developing  the  Model  II 


ft 
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correction  [Hopfield,  1965].  A  two-part  exponential  treatment  of  the  profile 
had  already  been  suggested  by  B.  R.  Bean  [1961],  with  especial  reference 
to  the  bending  of  a  low -angle  radio  ray  in  tropospheric  propagation.  As 
regards  satellite  use,  however,  the  bi -exponential  has  the  same  disadvantages 
as  a  simple  exponential  profile  and  was  therefore  not  adopted.  The  recent 
three -parameter  exponential  treatment  [Bean  et  al,  1966]  could  provide  the 
proper  area  under  the  N  profile;  but  with  added  problems  in  evaluating 
the  correction.  It  has  not  been  used  here. 

It  was  found  practical,  for  satellite  doppler  or  range  data,  to 
employ  a  two-part  treatment  which  uses  two  functions  of  the  form  of  equation 
(5),  with  p,  =  4  for  each  part.  The  resulting  "dry"  profile  is 

Nd  *  “d  \  - h)k  ■  hs  V 

a  a 

If  h  is  chosen  to  give  the  observed  area  under  a  "dry"  profile  and  the 
°d 

value  of  kj  is  determined  from  this  and  the  surface  conditions,  the  re¬ 
sulting  profile  gives  a  good  approximation  to  the  shape  of  the  obser¬ 
ved  profile  up  to  the  usual  height  of  balloon  data,  even  the  lapse  rate 
near  the  surface  ordinarily  being  within  a  few  percent  of  the  observed 

value.  The  value  of  h  is  of  the  order  of  1*0  km  and  will  be  discussed 

°d 

below. 

A  quartic  equation  is  also  a  reasonably  good  approximation  to  the 
usual  decrease  of  Nw  with  height,  i.e.. 
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N  =  k  (h  -  h) 
w  w  '  o  . 

w 


h  s:  h 


■w 


the  value  of  k  being  likewise  obtained  from  surface  conditions  and  a 
w 

known  or  postulated  h  ;  the  latter  ie  of  the  order  of  12  km  and  will 

°w 

also  be  discussed  below. 

The  total  theoretical  N  profile  is  then  the  sum  of  these  two 
components,  and  the  complete  expression  becomes: 


N 


N 


i 


E  * 


V 


¥ 


\ 


N 


i 


0 


h  >  h 

o 


J 


where  i  ■  1,2  refers  to  the  "dry"  and  "wet"  components  respectively.  All 
heights  are  measured  above  the  geoid.  The  subscript  T  refers  to  the 
tracking  station. 

Samples  of  observed  N  profile  components  (dry  and  wet),  and 
quartic  profile  components  designed  to  match  them  in  total  area,  are  shown 
in  Figures  l(a),  2(a),  3(a)  and  4(a)  for  quite  different  latitudes  and 
climates.  They  are  arranged  in  order  of  latitude,  heading  northward. 

The  total  profiles  are  given  in  the  (b)  curves  of  these  figures,  the  CRFL 
exponential  profile  being  also  shown  in  each  case  for  comparison.  The  shape 
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W)  TOTAL  N  PROFILE 

—  THEORETICAL  TWO-QUARTIC  PROFILE 

—  CRPL  EXPONENTIAL  PROFILE 
•  OBSERVED  DATA 


H*  1 


HEIGHT  (km) 

Fig.  1  PROFILE  OF  TROPOSPHERIC  REFRACTIVITY  N,  CANTON 
ISLAND,  LATITUDE  3°  SOUTH  >  JANUARY  1964 
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Fig.  2  PROFILE  OF  TROPOSPHERIC  REPRACTIVITY  N, 

WEATHER  SHIP  E,  LATITUDE  35°  NORTH -JULY  1967 
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Fig.  3  PROFILE  OF  TROPOSPHERIC  REFRACTIVITY  N, 

WASHINGTON,  O.C.,  LATITUDE  39*  NORTH  -  MAY  1964 
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of  the  two-quartic  expression  is  a  good  approximation  to  the  observed  data 
in  all  cases.  The  area  under  the  exponential  profile  is  too  small  by  only 
3  percent  for  the  Anchorage,  Alaska  profile  (Figure  4),  but  by  17  percent  for 
the  moist  tropical  climate  of  Canton  Island  (Figure  l). 

Observed  profiles  are  regular  in  form  and  differ  from  each 
other  to  a  relatively  minor  extent,  while  observed  Nw  profiles  may  be 
irregular  and  differ  widely  in  different  localities  and  seasons,  or  even 
from  day  to  day.  This  variation  will  be  discussed  in  a  later  section. 

IV.  Height  Parameters 

The  heights  h  and  h  of  Figures  1,  2,  3  and  4,  sel  cted  in- 
°d  °w 

dividually  for  profile  area-matching,  are  not  the  same.  Before  the  new  N 

model  can  be  used  routinely  to  compute  a  doppler  correction,  information  is 

needed  on  height  parameters  for  both  "dry"  and  "wet"  terms  for  any  location. 

Balloon  data  show  that  the  height  of  any  constant  pressure  surface 

in  the  atmosphere  is  to  seme  extent  a  function  of  latitude.  Hie  effective 

height  (for  refraction)  of  the  uncharged  part  of  the  atmosphere  is  probably 

also  a  function  of  latitude,  and  an  equatorial  bulge  is  likely. 

Initially,  latitude  variations  of  both  h  and  h  were 

°d  °w 

postulated.  An  attempt  was  made  to  find  values  of  h  and  h  at  the 

°d  °w 

equator,  and  the  latitude  variation  of  each,  using  a  least  squares  procedure 
to  match  theoretical  to  observed  values  of  jNdh  .  The  initial  results 
were  questionable,  and  a  different  approach  was  tried,  as  follows. 

The  simpler  assumption  was  made  as  a  first  approximation,  that 
the  height  hQ  is  invariant  with  latitude;  but  that  hQ  at  a  tracking 
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station  latitude  qp^,  is  given  by 


h 

°d(eq) 


+  A,  sin 


<&r 


(8) 


where  h  is  the  "dry"  height  at  the  equator  and  A,  the  amplitude  of 

°4(eq)  4 

variation  of  h  with  latitude.  If  h  has  an  equatorial  bulge.  A. 

°d  °d  4 

must  of  course  be  negative.  Integrating  equations  (7)  through  the  tropo¬ 
sphere,  the  total  area  under  the  two-quart ic  profile  is 


hQ 

J  dNdh  =  —  [N_  (h  -hJ+N  (h  -  hj]  (9) 

o  J  d  d  w  w 


Using  the  value  of  h  from  equation  (8),  and  a  postulated  value  of  h  , 

°d  °w 

in  equation  (9),  and  equating  the  resulting  theoretical  jNdh  to  a 

corresponding  observed  J*Ndh  from  balloon  data  (as  described  above),  the 

unknown  parameters  h  and  A.  were  obtained  by  a  least  squares 

°d(*q)  4 

procedure.  The  computation  was  repeated  several  times  with  different 

postulated  values  of  h  .  The  data  set  included  widely  varying  climates 

°w 

(e.g.,  Alaska,  Algeria,  Antarctica),  and  station  heights  from  sea  level  to 
1500  meters. 

The  resulting  parameters  are  tabulated  below  for  a  few  values  of 

h  ,  along  with  the  polar  values  of  h  which  result  from  these  parameters. 
°w  °d 
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TABLE  I 


Height  Pari  to  ter  s  for  Two-Quart  ic 

N  Profile  (tan) 

h 

h 

Polar  value  of 

°w 

°d(eq) 

a 

10 

43.858 

-5.986 

37.872 

12 

43.130 

-5.206 

37.924 

14 

42.402 

-4.426 

37.976 

Interestingly,  the 

polar  value  of  h 

is  nearly  the  same  i 

any  of  these  values  of  h  .  The  moisture  content  of  the  a  Jr  is  so  small 

°w 

in  polar  regions  that  the  solution  for  the  polar  h  is  not  much  affected 

°d 

by  assumptions  about  h 

°w 

Ihe  residual  errors  in  J*Ndh  (theoretical  minus  observed  value) 
have  an  rms  value  of  approximately  2  percent  for  any  of  these  sets  of  para¬ 
meters,  i.e.,  for  h  ranging  from  10  to  14  or  even  to  16  tar..  Thus  the 

°w 

area  (  ^Ndh)  alone  is  not  a  sufficient  criterion  for  deciding  between  the 

several  sets  of  parameters.  Comparison  of  the  shape  of  a  variety  of  observed 

with  theoretical  N  profiles  (from  equations  (7))  indicates  that  a  value 

of  h  =12  tan,  though  not  correct  in  all  cases  (e.g..  Figure  2),  is 
°w 

generally  satisfactory  and  is  suitable  for  preliminary  use.  From  Table  I, 
we  then  have 


h  =  43.130  -  5.206  sin*  <pT 

°d  i 


tan. 


(10) 
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at  a  latitude  <p>p  . 

The  ratio  of  theoretical  to  observed  jNdh  for  this  model  is 
plotted  in  Figure  5  as  a  function  of  total  surface  refractivity.  The  new 
N  model  makes  this  ratio  essentially  unity  for  any  value  of  surface  refrac¬ 
tivity,  and  unlike  Model  II,  needs  no  correction  factor. 

The  same  ratio  is  plotted  against  station  latitude  in  Figure  6. 
Any  systematic  departure  from  a  ratio  of  unity  is  small  enough  to  be  in 
doubt.  It  is  possible,  however,  that  slightly  different  height  parameters 
should  be  used  for  the  northern  and  southern  hemispheres  (in  analogy  to  the 
asymmetry  of  the  earth's  gravity  field).  No  further  work  has  been  done  on 
this  as  yet. 

Some  other  possible  systematic  errors  will  be  mentioned  later. 


V.  Tropospheric  Contribution  to  Satellite  Doppler  or  Range  Data 
A.  Theoretical  Expressions 

It  follows  from  equations  (l)  to  (3)  and  the  assumptions  about 
the  troposphere  that  the  tropospheric  contribution  to  the  doppler  shift  of 
a  satellite  signal  at  any  instant  is 


Af. 


tro 


3  (Ap 


tro 


SE 


(11) 


where  E  is  the  elevation  angle  of  the  slant  range  vector  to  the  satellite 
at  the  tracking  station  and  E  its  time  derivative  [Hopfield,  1963]. 

Using  the  procedure  developed  for  the  Model  I  correction,  but  the 
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N  profile  of  equations  ('j) ,  the  tropospheric  doppler  effect  becomes 


"tro*  10 


-6  |  rT  i  £  (\  \  <*> 

i-1,2 


(12) 


'i  i 


where  the  function  F.  (E)  for  either  component  (subscripted  "4"  to  refer 

4i 

to  the  quart ic  N  profile)  is 


F,  (E)  «  cos  E 
4i 


tro. 


*3  .  t  3  3  r2 

(i2  +-T 


tro 


i) 


2  3rTrtroi 

*1  (12  -  — 2— 1  +  3  r*tr0  > 


M13) 


3  rt~.  V  , 

+  < - T—  +  r\ro 


In  this  expression,  rT  and  rtrQ  are  distances  from  the  center  of  the 
earth  respectively  to  the  tracking  station  and  to  the  "top"  of  the  troposphere 
(dry  or  wet  component).  Also 


htroi  "  ho1  " 


rtro1  rT  +  htroi 


\  =  rT  sin  E 


=  rT  cos  E 


\  ^  -  */>1/2 
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The  value  of  F.  (E)  is  obviously  zero  at  E  =  90°  and  unity 
i 

at  E  =  0°  for  any  station  height  and  troposphere  height,  for  both  "dry" 
and  "wet"  terms. 

If  the  tropospheric  range  contribution  rather  than  a  doppler 
contribution  is  desired  (as  in  using  radar  or  integrated  doppler  data), 
this  is  obtained  by  evaluating  J(n-l)dp  with  the  same  input  information 
which  has  Just  been  described.  Hie  range  effect  (assuming  no  ray  bending) 
is  given  by  the  sum  of  the  "dry"  and  "wet"  contributions: 


Ap 


tro 


Y.  APi 

i-1,2 


where 


15  ll  -l  rT  rtro.  *1U12+  I 


,  _  2,31  3  o’io  s 

tr°i  1  2  troi  rT  '  3  Ttro1  i3±  +  15 


+  |  r.  2  (X  3  L  Jtp2)-  r  2  j (i  2  -ir  2) 
4  tro.,  3^  2  3,  2  tro.^  3±  3±  2  tro.. 


'i  3i 


♦  i  Ttro_,  ^ ♦  O  *  ^ 


troi  "3i  ) 

In  practice,  using  an  IBM  7094  computer,  the  computation  of  F.  (E) 

i 

or  <_\pi  is  done  in  double  precision  to  avoid  excessive  rounding  errors.  To 
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save  computing  time,  a  table  of  F.  (E)  vs  E  can  be  computed  for  each 

1 

satellite  pass,  since  r,p  and  r^rQ  pertain  to  a  particular  station;  the 
value  of  F>  (E)  for  a  specified  elevation  angle  can  then  be  obtained  by 
interpolating  in  the  table. 

If  a  tropospheric  range  rate  effect  is  desired,  this  is  obtained 
directly  from  equations  (12)  and  (13)  by  deleting  the  scaling  factor  — 
in  equation  (12). 

Equations  (12),  (13)  and  (l4)  come  from  exact  integral  expressions. 
They  are  based  on  physical  approximations  but  contain  no  mathematical  approx¬ 
imations.  To  the  extent  of  validity  of  the  physical  assumptions,  they  are 
therefore  good  at  all  elevation  angles,  and  do  in  fact  remain  finite  at  all 
angles.  The  tropospheric  contribution  to  range  rate  for  a  moderately  high 
pass,  obtained  from  equation  (l 2),  is  of  the  order  of  2  meters/sec  at  the 
geometrical  horizon,  where  the  actual  range  rate  is  approximately  6  km/sec, 
for  satellites  1000  km  above  the  earth.  The  range  contribution  from 
equations  (l4)  is  nearly  100  meters  at  the  horizon. 

The  effects  given  by  the  above  expressions  are  expected  to  be 
somewhat  larger  than  the  observed  ones  near  the  horizon  because  signal  path 
curvature  was  neglected.  The  actual  curved  path,  though  geometrically 
slightly  longer  than  the  slant  range  vector,  is  "optically"  shorter  (cf. 
Fermat ' s  principle ) . 

An  estimate  of  the  curvature  effects  neglected  in  equations  (12) 
and  (l4)  has  been  obtained  by  tracing  rays  through  sample  atmospheres 
described  by  equations  (7).  Preliminary  results  show  that  neglect  of  path 
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bending  in  a  nominal  atmosphere  may  make  the  computed  tropospheric  contri¬ 
bution  to  the  doppler  shift  or  range  rate  15  or  even  20  percent  too  large 
at  the  geometrical  horizon;  but  less  than  3  percent  too  large  at  elevation 
angles  above  5°.  The  neglected  path  curvature  error  in  the  range  effect  is 
approximately  half  as  large  as  this.  These  effects  are  shown  in  Figure  7 
for  humid  summer  conditions;  a  fairly  extreme  case.  More  details  will  be 
published  in  a  later  paper.  The  presence  of  the  curvature  effect  in  low- 
angle  doppler  data  is  readily  observable,  as  will  be  shown  below. 

B.  Observed  Results 

Tropospheric  effects  can  be  observed  (l)  in  the  doppler  residuals 

(observed  minus  theoretical  doppler  shift)  obtained  when  no  troposphe-  Le 

correction  has  been  applied  to  the  observed  data;  and  (2)  in  station 

position  errors.  Both  types  of  observed  effects  will  be  shown  here. 

In  processing  satellite  doppler  data,  the  thoretical  tropospheric 

correction  is  computed  from  equations  (12)  and  (13)  at  each  data  point  and 

applied  to  the  data  before  final  computations  are  made.  Inputs  are  geometric 

conditions  obtained  from  a  preliminary  orbit,  station  latitude  for  computing 

tropospheric  height,  and  local  surface  weather  data  observed  near  the  time 

of  the  satellite  pass.  Should  weather  data  be  lacking,  seasonal  values  of 

Nt  and  Nt  for  the  bracking  station  locality  are  used  ["Monthly  Climatic 
a  w 

Data  for  the  World;1'  also  Bean  et  al,  i960,  1966]. 

Figures  8  to  11  show  the  effectiveness  of  the  correction  described 
above.  Figures  8,  9  and  10  show  the  doppler  residuals  (observed  minus 
theoretical  doppler  shift)  for  three  satellite  passes  observed  under  differ¬ 
ent  geometrical  and  climatic  conditions:  in  each  case  (a)  without  and 
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Fig.  11  NAVIGATION  ERROR  IN  STATION-TO-ORBIT  RANGE, 
SATELLITE  1967  34A,  1  AND  2  SEPTEMBER  1967 
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(b)  with  the  use  of  the  tropospheric  correction.  The  thoretical  correction 
is  also  shown  for  each  pass.  In  each  figure  the  uncorrected  residuals  of  part 
(a)  follow  the  trend  of  the  thoretical  tropospheric  effect.  Little  sys¬ 
tematic  error  remained  after  the  tropospheric  correction  had  been  applied 
(part  (b)  of  each  figure).  In  the  one  case  where  actual  tracking  extended 
below  the  geometrical  horizon  (Figure  10),  the  magnitude  of  the  theoretical 
correction  was  noticeably  too  large  below  2°  elevation.  This  is  compatible 
with  the  expected  but  disregarded  path  curvature  effect  shown  in  Figure  7- 
Hie  effect,  however,  is  a  function  of  local  conditions  and  its  magnitude 
would  often  be  less  than  that  of  Figure  7. 

It  should  be  noted  that  Figures  8  to  10  show  doppler  residuals  for 
passes  in  which  non-tropospheric  errors  were  very  small.  Other  errors  can 
easily  mask  tropospheric  effects  in  the  central  portion  of  a  pass,  but 
even  then,  the  characteristic  large  increase  of  the  tropospheric  effect 
near  the  horizon  is  routinely  visible. 

The  effect  of  the  troposphere  on  station-to-orbit  range  is  shown 
in  Figure  11.  (Navigated  station  position  in  a  direction  parallel  to 
the  orbit  is  affected  by  the  troposphere  only  if  the  tracking  data  are  not 
symmetrical  about  the  point  of  closest  approach;  and  will  not  be  discussed 
here.)  Figure  11  shows,  for  passes  at  different  elevation  angles,  the 
errors  in  station-to-orbit  ranges  ("navigated"  ranges  as  compared  with  the 
best  available  theoretical  values)  for  a  two-day  set  of  satellite  passes, 
observed  at  tracking  stations  distributed  over  the  earth.  Each  point  repre¬ 
sents  one  pass.  Parts  (a)  and  (b)  show  the  range  errors  obtained  from  the 
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same  tracking  data,  respectively  without  and  with  the  use  of  the  tropospheric 
correction;  using  all  tracking  data  obtained  down  to  an  elevation  angle  of 
1°. 

Tropospheric  effects,  when  uncorrected  as  in  part  (a)  of  figure  11, 
shorten  the  apparent  station-to-orbit  range  by  20  meters  (average)  in  the 
case  of  high-angle  passes,  but  by  very  much  more  for  lower  posses  (40  meters 
for  a  pass  whose  maximum  elevation  is  30°,  100  meters  for  a  10°  pasB,  and 
more  than  that  for  still  lower  passes). 

Hiese  systematic  errors  were  removed  by  the  use  of  the  tropospheric 
correction,  as  shown  in  part  (b).  For  passes  lower  than  10°,  however,  the 
average  correction  was  a  little  too  large.  This  over-correction  can 
reasonably  be  ascribed  to  the  (neglected)  low-angle  path  curvature  effect 
which  has  already  been  pointed  out  in  the  doppler  residuals. 

The  over-correction  disappeared  when  data  were  deleted  below  an 
elevation  of  5°,  as  in  Figure  ll(c).  Here  no  appreciable  systematic  effect 
remains.  The  correction  was  satisfactory  on  the  average,  and  the  neglected 
effect  of  signal  path  curvature  was  apparently  unimportant  as  long  as  the 
satellite  elevation  during  the  data  interval  was  at  least  5°* 

VI.  Discussion 

The  residual  range  errors  of  Figure  ll(c)  show  little  or  no  bias 
as  a  group.  Their  scatter  has  many  possible  sources,  including  uncorrected 
tropospheric  effects.  Such  effects  may  produce  an  error  in  the  range 
determination  which  is  systematic  for  any  one  satellite  pass.  lhe  fact 
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that  the  scatter  is  noticeably  greater  for  the  lower-angle  passes  indicates 
the  atmosp'.iere  as  one  probable  source. 

On  the  basis  of  probable  errors  in  the  inputs  and  assumptions  used 
above,  it  is  estimated  that  for  any  one  pass,  tropospheric  effects  in  doppler 
data  are  reduced  at  least  by  an  order  of  magnitude  by  the  correction  given 
above.  A  few  remarks  will  now  be  made  about  the  remaining  errors,  and  about 
some  areas  for  further  investigation;  in  particular,  height  parameters,  and 
local  weather  effects. 

The  area  under  any  Ny  profile,  even  for  a  humid  atmosphere, 
represents  a  small  fraction  (seldom  more  than  10  percent)  of  the  area 
under  the  total  N  profile  (Figures  1  to  4);  and  contributes  only  this 
same  fraction  to  the  total  tropospheric  effect  (range  or  doppler),  in  the 
zenith  direction.  Hie  same  rule  holds  approximately,  down  to  about  5° 
t  ''evation,  but  not  below.  At  the  geometrical  horizon,  the  moisture  content 
of  the  air  may  account  for  as  much  as  25  or  30  percent  of  the  total  tropospheric 
doppler  effect  (from  equation  (l 2)).  If  data  below  5°  elevation  are  deleted, 
however,  the  "dry"  portion  of  N  accounts  for  90  percent  or  more  of  the 
tropospheric  range  or  doppler  effect. 

Thus  an  error  of  20  percent  in  representing  the  Nw  profile  area 
changes  the  total  tropospheric  effect  by  2  percent  or  less,  while  an  error 
of  only  2  percent  in  the  profile  area  changes  the  total  effect  by 

nearly  2  percent.  The  relatively  minor  importance  of  the  "wet"  term,  for 
satellite  range  or  doppler  applications,  which  was  not  recognized  when  the 
Model  I  troposphere  correction  was  developed,  became  obvious  when  the  two 
profile  components  were  studied  separately.  For  future  improvement,  the 
height  parameters,  especially  for  the  "dry"  term,  should  be  better  determined 
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as  to  latitude  variation,  local  climatic  characteristics,  and  possibly  time 
variation  (e.g.,  seasonal  or  diurnal). 

An  illustration  of  local  weather  effects  is  presented  in  Figures 
12  and  13:  two  sample  N  profiles  obtained  from  individual  balloon  ascents 
at  Pensacola,  Florida,  six  days  apart;  with  their  companion  theoretical  pro¬ 
files.  The  theoretical  "dry"  height  was  routinely  obtained  from  equation  (10), 
the  assumed  "wet"  height  being  12  km.  The  observed  "dry"  component  profiles 
are  relatively  smooth  and  similar  in  shape,  but  there  is  a  marked  difference 
between  the  "wet"  components  for  the  two  days.  These  "wet"  profiles  can  be 
qualitatively  correlated  with  weather  maps  for  the  dates.  A  cold  front 
had  reached  the  location  shortly  before  the  Figure  12  observations.  The 
arriving  air  mass  was  "continental  polar";  however,  this  had  not  yet  reached 
very  high,  as  can  be  see.  the  increase  of  Ny  with  height,  just  uoove 
the  surface.  At  the  time  of  the  Figure  13  observations,  there  had  been 
no  recent  weather  front  and  the  air  mass  in  the  region  was  "maritime  tropical" 
and  showed  a  steep  decrease  of  Nw  with  height.  Minor  irregularities  in  the 
"wet"  terms  occurred  at  higher  altitudes  on  both  days.  In  spite  of  the  great 
difference  in  the  "wet"  profiles,  the  observed  value  of  the  zenith  integral 
jNdh  (total)  was  only  one  percent  higher  for  Figure  13  than  for  Figure  12. 

A  two-quartic  correction  based  on  surface  data  would  have  been  3  percent 
higher  for  the  Figure  13  than  for  the  Figure  12  data,  leaving  a  2  percent 
differential  in  the  precision  of  the  correction  for  the  two  days. 

Aside  from  temporary  weather  variations  and  systematic  latitude 
variations,  other  climatic  differences  may  exist  for  which  some  systematic 
allowance  could  be  made.  Temperature  Inversions,  for  example,  are  so 
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Fig.  12  REFRACTIVITY  PROFILE,  PENSACOLA,  FLORIDA 
21  APRIL  1965,  12  HOURS  U.T. 
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common  as  to  be  expected  in  some  areas  and  seasons.  If  the  optimum  correc¬ 
tion  is  desired  at  a  specific  tracking  station,  the  N  profile  in  the 
region  should  be  sufficiently  observed  to  determine  its  systematic  and 
predictable  variations  from  the  nominal  profile. 

The  N  profile  presented  above  may  be  of  interest  also  for 
correcting  optical  data  [Guier,  1968],  since  the  "dry"  term  of  the  radio 
refractivity  is  very  nearly  the  same  as  the  optical  refractivity  of  air. 

The  quartic  "dry"  component  profile  could  be  used  with  little  modification 
for  laser  and  other  optical  applications. 
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